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We develop a strong coupling approach for a general lattice problem. We argue that this strong 
coupling perspective represents the natural framework for a generalization of the dynamical mean 
field theory (DMFT). The main result of this analysis is twofold: 1) It provides the tools for a unified 
treatment of any non-local contribution to the Hamiltonian. Within our scheme, non-local terms 
such as hopping terms, spin-spin interactions, or non-local Coulomb interactions are treated on equal 
footing. 2) By performing a detailed strong-coupling analysis of a generalized lattice problem, we 
establish the basis for possible clean and systematic extensions beyond DMFT. To this end, we study 
the problem using three different perspectives. First, we develop a generalized expansion around the 
atomic limit in terms of the coupling constants for the non-local contributions to the Hamiltonian. 
By analyzing the diagrammatics associated with this expansion, we establish the equations for a 
generalized dynamical mean-field theory (G-DMFT). Second, we formulate the theory in terms of 
a generalized strong coupling version of the Baym-Kadanoff functional. Third, following Pairault, 
Senechal, and Tremblay— , we present our scheme in the language of a perturbation theory for 
canonical fcrmionic and bosonic fields and we establish the interpretation of various strong coupling 
quantities within a standard perturbative picture. 

PACS numbers: 71.10.Fd 



Understanding strongly correlated electron systems represents a major challenge in solid state physics. Dynamical 
mean-field theory (for a review see RefQ) has emerged as a powerful tool to address this challenge. Within this 
approach, the lattice problem for an interacting electron system is mapped onto a quantum impurity problem with 
a "bath" determined self-consistently^. In its simple form, the construction becomes exact in the limit of infinite 
lattice coordination as shown in Ref. U In addition, several generalizations of the approach, such as the extended 
DMFT 4 ! 5 ! 6 ! 7 ! 8 and the DMFT treatment of correlated hopping^ifl, requiring different scalings of the parameters in 
the Hamiltonian, have been presented. 

More generally, DMFT is a tool with strong computational power, and provides valuable insight into the physics 
of strongly-correlated electrons. In the case of simple model Hamiltonians with purely local interactions, such as 
the Hubbard model, it was able to provide a consistent description for several non-perturbative properties such as, 
for example, the Mott-Hubbard metal insulator transition!. Currently, many efforts are made to extend the DMFT 
approach in order to make it suitable in addressing more realistic problems. 

The purpose of this work is to develop a strong coupling approach for interacting lattice problems. Within this 
approach, we construct a systematic generalization of the DMFT technique capable of describing any non-local term 
of a generic lattice Hamiltonian. The present formulation has several advantages. First, it unifies various DMFT 
schemes. For example, the treatment of correlated hopping2ii&, the single site DMFT of the Hubbard model (see 
Section III), and the extended DMFTiiSffl are all limiting cases of this unified approach. Second, besides unifying 
different DMFT schemes, this formulation is sufficiently general to be a starting point for expansions around DMFT. 
One possibility is described in the Appendix. In addition, the formulation in terms of Hubbard operators makes it an 
ideal framework for a DMFT based renormalization group approach along the lines of Ref. ITil Third, being based 
on an expansion around the atomic limit, it provides new and valuable insight into the diagrammatics underlying 
cluster DMFT schemes, which will be a subject of a subsequent publication. 

The central result of this work is the formulation of a perturbation theory, for a general Hamiltonian of the form 



containing a local term, Hq, and a non-local one, H±. In terms of Hubbard operators, X. , we can write the two 
contributions as 



I. INTRODUCTION 



H — H + Hi , 



(1) 




(2) 



a: 
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and 

H! = J2 J2 E^'xfxfr (3) 



In these equations a, j3, a' and j3' represent single-site states. For simplicity, we consider the case of spin | fermions, 
so that for each site i there will be four states \a): |0) (empty site), | f) (single occupied site with spin up), | J.) (single 
occupied site, spin down), and 2) (double occupied site). The Hubbard operators X^ describe transitions between 
these states. These operators have a fermionic character if the occupation numbers of the states a and (3 differ by 
one, and bosonic character otherwise. The algebra of the Hubbard operators is defined by the multiplication rule 

Xfxf a '= 6^xr\ (4) 



together with the conserving condition 



and the commutation relation 



E*r = i: (5) 



[Xfxf a '}± = S^iS^xr' ± 8 aa .xf p ). (6) 

In Eq. © the anticommutator is used if both operators are fermionic and the commutator otherwise. The canonic 
fcrmion operators Ci a can be expressed in terms of Hubbard operators as 

c la = X?° + aXf 2 . (7) 

The parameters A Q in Eq. © represent the single-site energies and will be determined, in general, by the chemical 

potential, /x, the on-site Hubbard interaction, U , and the external fields. The coupling constants E"^ a 13 may include 
contributions coming, for example, from hopping, £y, spin-spin interaction, Jy, or non-local Coulomb interaction, 

Vij. Within the present approach, all th 
generalized "hopping" matrix elements. 



Vij. Within the present approach, all these contributions are treated on equal footing, and we can regard E°f a ^ 



Our main object of interest is the Green's function of the Hubbard operators. However, the theory is expressed 
naturally in terms of a generalized irreducible two-point cumulant, M, and a dressed hopping, £. Using the language 
of a standard field theory (see the Appendix) , we interpret M as a generalized self energy and £ as the corresponding 
Green's function associated with a set of auxiliary canonical fermionic and bosonic fields, and show that they have a 
clear diagrammatic interpretation in the locator expansion. In addition, our formulation contains the average Q = (X) 
of the Hubbard operator, which can be viewed as a generalized "magnetization" , together with a generalized effective 
"magnetic" field, h, represented, in the language of the standard field theory, by the mean value of the auxiliary field. 
Formally, the relationship between these quantities and the central object of the theory, G, is given by the equations: 

G = ((XX)) - (X)(X) = (M- 1 - E)-\ 

and 

£ = E{1- ME)- 1 . 

where E represents the bare coupling constants. 
Next, we introduce the functional 

T[£, Q] = -Tr ln[l - ME] - Tr[M£] - ^QEQ + Q], 

where ^ represents a generalized Baym-Kadanoff-type functional that can be obtained as a sum of all vacuum-to- 
vacuum skeleton diagrams. Using this functional perspective, we can view £ and M, as well as Q and h, as pairs of 
conjugate quantities and we have 6 1 $>/8£ = M and 6"i>/SQ = h, where the effective "magnetic" field can be expressed 
as h — EQ. 

The single site dynamical mean field theory is a local approximation for M w M; oc , which can be obtained as a 
stationarity condition for the functional T. This condition translates into a self-consistent impurity problem defined 
by the statistical operator 

ft mp = e-P H o mp f exp J - J X{t)A(t - t')X(t') + h imp {T)X( T ) 1 , 
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where hi mp represents an external field and A the hybridization. These quantities are subjected to the self-consistency 
conditions: 

^ E KK^n) - E{k)]- 1 = [M," V«) - A(iu n )}-\ 

k 

and 

h imp + A(0)(X) =E(X). 

The strategy that we use in presenting our G-DMFT scheme consists in constructing three formally distinct formu- 
lations. On the one hand, this allows us to establish the equivalence of these approaches to the problem of correlated 
electrons and to extract a generalized unitary picture. On the other hand, these various angles clearly reveal the 
aspects that are relevant for the possible extensions of the theory. The natural starting point for our construction 
is given by a generalized expansion around the atomic limit. This expansion, in terms of the coupling constants 
of the non-local contributions to the Hamiltonian, is derived in Sec. II. The technique is characterized by a close 
formal analogy with the "canonical" perturbation expansion in terms of inter-site hopping. Consequently, we will 
focus on the new aspects generated by the description of the problem in terms of Hubbard operators, the main goal 
being to obtain the relation between the Green's function for the Hubbard operators and the irreducible two-point 
cumulant. Using this relation, we derive in Sec. Ill the generalized DMFT equations. A simple illustration of the 
implementation of our scheme is given for the Hubbard model. In Sec. IV we formulate the theory in terms of a 
functional of the renormalized "hopping" ' The generalized DMFT equations are shown to be the result of a simple 
local approximation on a generalized Baym-Kadanoff-type functional. An alternative derivation is described in the 
Appendix. By decoupling the non-local term via a Hubbard-Stratonovich transformation, we formulate the theory in 
terms of a set of canonical fermionic and bosonic fields. In this language, the standard treatment based on Wick's 
theorem applies. A possibility of expanding around the DMFT solution is also presented. 



II. EXPANSION AROUND THE ATOMIC LIMIT 



Expansions around the atomic limit have been powerful tools for studying both models with localized spins 15 , 
and models of interacting itinerant fermions. Following the pioneering work of HubbarcU&, Metzne»i£ developed 
a renormalized series expansion for the single-band Hubbard model. Generalizing the "linked-cluster expansion" 
ideasA^, this approach involves only connected diagrams and unrestricted lattice sums, and enables one to construct 
self-consistent approximations. By analyzing this renormalized expansion, it was showni that, in the limit of infinite 
spatial dimensions, the dynamical mean-field theory equations for the Hubbard model are recovered. An extension 
of this approach to the case of correlated hopping was developed by Shvaika 10 . The purpose of this section is to 
obtain a generalization of the DMFT equations that describe the physics of a generic lattice model for interacting 
fermions, starting from a renormalized strong-coupling expansion around the atomic limit 11 . Within this scheme, all 
the non-local contributions to the Hamiltonian, for example hopping terms, spin-spin interaction terms, or non-local 
Coulomb interaction contributions, are treated on equal footing. 

Let us consider a system described by the Hamiltonian The first step in our derivation is to write an expansion 
for the grand-canonical potential and the Green's functions for the X-operators in terms of "hopping" matrix elements 
E"^ a ^ and bare cumulants. Starting with the statistical operator 

p = e-P Ho a{p) =e-P Ho f exp j- jf" drtfi(r) j , (8) 
we can write the grand-canonical potential as 

n = n -iln<a(/3)) , (9) 

where £Iq — —4 \nTr(e~^ H °) is the grand-canonical potential in the atomic limit, (3 ~ l/ksT is the inverse tempera- 
ture, and the unperturbed ensemble average is given by 



( . . . } = Tr{exp(-0H o ) . . . }/Tr{exp(-/3ff )}. 



(10) 
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Expanding the exponential in Eq. JSJ) we obtain the nth order contribution to (a(/3))o: 



(HP))o l) = I' dn... I dr^fiHxin) . . . H n (r n )}) 



P rP 



P 



^ ^ EgfroiPt . . . Et:f n na ^ / dn...dr n (11) 

{ik,jk}{a k ,Pk,a' k ,P' k } J ° 

x (f{X^(n)X^( n ) . . . X^(r n )X^(r n )}) Q , 

where T represents the imaginary-time ordering operator, and X" (r) are the Hubbard operators in the interaction 
representation, 

Xf{ T ) = e HoT xfe- HoT . (12) 

Using the multiplication properties Q of the Hubbard operators, it is straightforward to obtain the explicit imaginary 
time dependence, 

Xf(r) = e -^"- x ^ T Xf. (13) 

Because Hq is a sum of local operators, the ensemble average in Ijllf) factorizes into independent local averages that 
can be evaluated using the algebra of the X-operators. However, the explicit calculation of f2 is cumbersome, as the 
summation over the site variables it is restricted. To overcome this problem, following Metzner— , we introduce the 
bare cumulant defined by 



S n W 

@i( a ufhi ■ ■ ■ > a n,0n;n, ■ ■ ■ ,T n ) 



(14) 

{r,=0} 



5"na l pA T l) ■ ■ ■ dVa n P n (Tn) 

where the generating functional is 

W[{r ? }] = ln/fe X pJ -J2 dr rj^X^iA \ . (15) 



a/3 



The fields rj a /s are either complex numbers or Grassmann variables, depending on the nature (bosonic or fermionic) 
of the corresponding operators X a @. In the Eqs. 1)14(1 and (|15[1 the site index, i, for the r\ field and the Hubbard 
operators, was omitted for simplicity. 

Using the bare cumulants allows one to express the grand-canonical potential as an expansion containing unrestricted 
sums. Each term will be a product of "hopping" matrix elements, E^ a ^ , and local cumulants, and can be represented 
diagrammatically. Further, using linked-cluster type arguments^ one can show that only the connected diagrams 
contribute. The basic diagrammatic rules are similar to those given by Metzner— for the Hubbard model. However, 
the internal lines corresponding to "hopping" matrix elements will carry two extra pairs of indices (af3) and (a 1 (3') 
(one for each end of the line) representing the single-site states. Also, the nth order cumulant, C- \ represented 
by a n-valent point vertex, will be attached to n lines, with the corresponding labels {otkPk)- If we split the states 
associated with a cumulant in two subsets, (a%, . . . ,a n ) and (/3\, . . . ,[3 n ), the multiplication rule Q requires that 
each single-particle state (|0), |er), or |2)) be found the same number of times in each of the subsets. Also, in contrast 
with the simple case of an expansion in the hopping matrix, £y , in general, cumulants of order n with n odd may be 
non-zero. This is the case if Hi contains terms with bosonic X-operators having a non-zero average value. Such terms 
may occur, for example, if a non-local Coulomb interaction VijUiUj is considered. To summarize these observations, 
we present in Fig. Q]the diagrams yielding the leading contributions to il. 

Our next task is to write an expansion for the Green's functions of the Hubbard operators. We define 

G a0 a 'p>(i,r,T,T') = -Tr{p f[Xf(r)Xf(r')]}, (16) 
where the statistical operator p is given by Eq. JSJ, and the Hubbard operators are in the Heisenberg representation, 

Xf{r) = e HT xfe~ HT . (17) 



The perturbation expansion for the Green's functions is similar to the expansion for the grand-canonical potential. 
The consequence of the presence of two extra X-operators in (|16fl is that the corresponding diagrams will be rooted, 
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FIG. 1: Diagrams contributing to the expansion of the grand canonical potential in terms of generalized "hopping" matrix 
elements (full lines) and bare cumulants (small circles). The contributions from diagrams a, b, d, e, f, and g are non-zero only 
if Hi contains terms with bosonic X-operators having non-zero average values (Q-operators). 

i.e. one (for i = j) or two of the vertices will have fixed site indices and some fixed single-site state indices. The 
general rules are analogous to those for the Hubbard model 17 , with the difference that extra indices (associated with 
the single-site states) are carried by the "hopping" matrix elements and the bare cumulants. A special case occurs 
when bosonic X-operators with non-zero average values, (X al3 ) = Q af3 , which we will call Q-operators, are present 
in Hi. In this case, as mentioned before, cumulants of order n with n odd may be non-zero, leading to disconnected 
rooted diagrams in the expansion of G a pa'P'- By simply regrouping the contributions from disconnected diagrams we 
can write: 

G a p a'fi'(i,r,T,T') = G aP a'/J'&ijr.O - (Xf){X? 0> ), (18) 

where {X"^ 3 ) represents the average of the operator X" 13 with respect to the full Hamiltonian, and G a p a 'P' bas a 
diagrammatic representation containing only connected diagrams. In the rest of this work we will always use the 
"connected" Green's functions defined by (JTHJ. Obviously, if Hi does not contain Q-operators, G and G coincide, as 
the averages in ()18|) vanish. To summarize, we present in Fig. [21 some leading contributions to G. 

At this point, a brief discussion about the Fourier transformation of the Green's functions and of the cumulants 
is required. Due to the fact that the Hubbard operators have either bosonic-like or fermionic-like characters, the 
Fourier transformed Green's functions, G a p a >pi (i, j; ioj n ), will have odd Matsubara frequency, if both X-operators 
are fermionic, or even Matsubara frequency otherwise. On the other hand, the cumulants may depend on both even 
and odd Matsubara frequencies. Recalling the definition, Eqs. I|14ll5p . a pair of indices otkftk associated with a 
fermionic-like Hubbard operator will induce a dependence on an odd frequency Wfc, while an even frequency will be 
associated with a pair with bosonic character. A line connecting two vertices, and yielding the factor Efj 3 a ^ , will 
carry the frequency corresponding to the nature of the two pairs (a/3) and (a' '/?'). This nature is always the same as 
a result of the fact that the Hamiltonian conserves the number of particles. In the rest of this work we will not use 
different notations for odd and even frequencies, but we should keep in mind this aspect of the problem. 

The expansion that we derived is rather complicated and it may be of little help if used directly. However, we 
know that in order to obtain self-consistent approximations and, in particular, DMFT-type equations, a renormalized 
expansion is necessary. The key in obtaining a renormalized expansion is to introduce the notion of irreducibility 

with respect to one line (representing Ef? a 13 ). We define the irreducible two-point cumulant, Mf? a 13 (io-i n ), as the 
sum of all connected diagrams having two external legs (labeled by a/3 and a'/3' , respectively), with the property that 
the two legs cannot be separated by cutting one line. This definition differs slightly from the standard defmitioniil, 
which requires that the diagram be not disconnected by cutting a single line. The difference comes from the existence 
of cumulants with odd number of legs. If Hi does not contain Q-operators, the standard diagrammatics, along with 
the usual definition of irreducibility are valid. Some typical, low order contributions to the irreducible cumulant 
MfP a 13 (iu> n ) are shown in Fig|3| Notice that we can classify the diagrams into classes, each class containing one 
"basic" diagram that is one-particle irreducible with respect to all "hopping" lines plus diagrams derived from it 
by decorating the bare cumulants with reducible contributions. It is useful to perform partial summations of the 
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FIG. 2: Diagrams representing the leading contributions to the expansion of the Green's function around the atomic limit. 
The cumulants characterized by fixed site and state labels have small legs attached to the corresponding circle. For n = 3 only 
diagrams containing cumulants with an even number of legs are shown. 
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FIG. 3: Diagrammatic representation of the renormalized expansion of the Green's function and some leading contributions 
to the irreducible two-point cumulant M??(w n )- The simplified notation (a/3) = fi, (a'/3') = f was used. The diagrams 
containing one-particle reducible decorations of bare cumulants can be summed and replaced by diagrams with decorations 
equal to ^ . v E^ - "(Xj). The average (Xj) is represented diagrammatically by a square. 



diagrams in each class. The resulting diagrams will have all the reducible decorations replaced by factors equal to 

Erpa/3 a'fj' / v a'/3'\ 

We can now write a renormalized expansion for the Green's functions in terms of irreducible cumulants and "hop- 
ping" matrix elements. The diagrammatic representation of this expansion is shown in FigO Analytically, the 
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corresponding equation will be: 

G% = Mr + M ^ E uGl»i (19) 

k,l P.7T 

where, to reduce the number of indices, we introduced a simplified notation for a pair of single-site states, (a/3) = fx. 
It is also convenient to regard a quantity A^ v as an element of a matrix A. Using this matrix notation, and Fourier 
transforming the spatial dependence, one can write the relation between the irreducible cumulant and the X-operator 
Green's function, 

G(k, iw n ) = [M- J (k, iw n ) - E(k)] _1 , (20) 

where E(k) is the Fourier transform of the coupling constant matrix. What exactly are the Green's function included 
in the matrix G? Obviously, those that are coupled by the generalized "hopping" £r. . Specifically, if {/i n }n=i,...,i 3 
and {v n }n=i,...,p are the sets of single-site pairs that label the coupling constants in Hi, each pair appearing only 
once in the corresponding set, the matrices will have dimension P and the matrix elements will be labeled G VnlJ " a , 
M" nfln , and E^ nVn , respectively. A particular example of this construction will be presented in the next section for 
the Hubbard model. 

Equations 1)19(1 and 120(1 represent the main result of this section. They are valid in arbitrary dimensions and for 
any Hamiltonian that has the form given by Eqs. ((113(1 when written in terms of Hubbard operators. 



III. GENERALIZED DMFT EQUATIONS 

In the standard dynamical mean-field theoryi, the single-particle self energy for the electron Green's function is 
local in the limit of infinite spatial dimensions. However, for a general lattice Hamiltonian containing fermionic X- 
operators, this property is violated, as shown by Schiller 9 in the case of correlated hopping. Nevertheless, one can 
overcome this problem by adopting a strong-coupling perspectiveAC* and regarding the two-point irreducible cumulant, 
rather than the self-energy, as the the quantity in which the theory should be formulated. If Hi contains only fermionic 
X-operators, and assuming that the coupling constants in Hi are properly scaled, to keep the energy (per site) finite, 
one can showiii that the irreducible cumulants become local, 

Mijfan) = S l3 M (iuj n ). (21) 

Consequently, the natural quantity in which such a theory should be formulated is the irreducible cumulant, rather 
than the self-energy. This becomes immediately transparent using Eqs. ((20(1 and 1(21(1 . and the fact that the electron 
Green's function is given by the sum 

Q a (k, iu; n ) = G 0a CT °(k, ioj n ) + G s2 2S (k, iu n ) + <j[G 0a 25 (k, ioj n ) + G s2 CT °(k, iw n )}. (22) 

Although Mo is k- independent, Q will have, in general, a complicated dependence on the wave-number, and the 
corresponding self-energy will be non-local, breaking the standard DMFT scheme. Formally, a self-energy matrix can 
be introduced by defining 

S(z) = (z + n)I - E - G-^z) = (z + m)I - M^z), (23) 

with I being the identity matrix. However, in general, this quantity does not have the desired analytic properties 10 
and may diverge as \z\ — > oo. 

If we go further and generalize the model to include both fermionic and bosonic degrees of freedom, the picture 
presented above fails, as one cannot scale the coupling constants in such a way that one has both finite energy (per 
site) and local irreducible cumulants in the limit of infinite dimensions^. In contrast, this work assumes a different 
perspective: rather than viewing DMFT as an approximation that becomes exact in the limit of infinite coordination, 
we see it as the first order approximation in a hierarchy of (increasing size) cluster approximations. In this picture, 
the exact limit is obtained for an infinite-large cluster. Consequently, our G-DMFT scheme will be defined by the 
condition that the two-point irreducible cumulant be local. Unlike the case of extended DMFT, we have no small 
parameter (other than the inverse of the cluster size) justifying this approach. Obviously, if the model contains only 
fermionic degrees of freedom, the two perspectives described above are equivalent . 

Based on the locality of M, we can approach the central step in a DMFT-type scheme, namely the self-consistent 
mapping of the lattice problem onto an effective impurity problem. Let us consider the single-impurity problem with 
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the statistical operator 



Pimp 



= e -PK mp fex V \-f dT f dr' E X afi {T)A aP a 'P\T-T t )X a 'P'{r')+ E h^ p {r)X afi (r) 

I /n a/3 a'p' a,/3£B Q 



(24) 

where hi mp (T) is an external field that couples to the Q-operators, Bq is the set of indices corresponding to the 
Q-operators, A(t — r') is the hybridization matrix and H l ™ p contains one term, for example i — 0, from the local 
Hamiltonian (01 . Using the expansion procedures described in the previous section, with the hybridization replacing 
the "hopping" matrix E, we obtain for the impurity Green's function a relation similar to (|2t)(l . 

G imp (iu n ) = [M~^ p (iuj n ) - A(iu) n )] , (25) 

where all the quantities are local, and M. imp is the irreducible cumulant with respect to the hybridization. For practical 
calculations, as well as to gain some physical insight, it is useful to express the impurity problem in a Hamiltonian 
form. Taking into account that in our matrix formulation the indices (uf3) corresponding to fermionic and bosonic 
Hubbard operators do not mix, all the matrices, including A Q/3 a @ , are block-diagonal with one block corresponding 
to the fermionic indices and one to the bosonic sector. Consequently, we will introduce two types of auxiliary degrees 
of freedom, one associated with a fermionic bath, described by operators (aka, a \ a ), an d one associated with a bosonic 
bath, described by (bka,b k(T ). The impurity Hamiltonian will be 

(a/3)eBq k,a k,a 

+ E E ^(k)al a( X^ + h.c.]+Y, E K (k)bl a ^ + h.c.}. (26) 

k (ul3)£F k (af3)eB 

The first term in the right hand side describes an isolated impurity, the second term represents the coupling to the 
external field and the summation is restricted to labels corresponding to bosonic operators having a non-zero mean 
value, the Q-operators (Bq), the third and forth terms describe the fermionic and bosonic baths, respectively, and the 
last two terms represent the coupling of the impurity with the two baths. In the last two terms, the summations over 
the X-operator labels (a/3) are restricted to the fermionic (T) and bosonic sectors (B), respectively. The expression 
of the hybridization in terms of coupling constants V a ^(k) can be obtained by integrating out the auxiliary degrees 
of freedom. We have 

A«* ^'(iu> n ) = A E V "" ma " {k \ (27) 

where I = a if the element is in the fermionic sector of the A matrix and I = b for the bosonic case. 

Returning to our main problem, we obtain the self-consistency conditions by requiring that the lattice irreducible 
cumulant be equal to the effective impurity problem irreducible cumulant, M imp (z) = M (z). By comparing the 
diagrammatic expansion of Mj„ lp in terms of bare cumulants (dressed with external field lines) and hybridization 
functions with the expansion of Mo, we observe that the equality is satisfied if we identify a line representing 
with the sum of all the paths starting with Eftf and ending with EJq , as well as the sum betweeen an external 
field line and the one-particle reducible decoration J2„ A A " y (X l/ ) with the lattice one-particle reducible decoration 
Y) j v E^- (Xj). Explicitly, the expansion for the impurity irreducible cumulant is identical with the expansion of the 
lattice irreducible cumulant provided 

A(iw„) = E E 0l G l ( ° ) (z^ n )E j0 , (28) 

where G^\iuj n ) represents the matrix Green's function for a lattice with the site "0" removed, together with a 
condition that determines self-consistently the external field h, 

h? mp = E E£(X»)-J2^(0)(XZ)- (29) 

j, v^Bq v 

Notice that equation (|28|l represents the generalization of the well known DMFT self-consistency condition 1 . In 
addition, if the non-local Hamiltonian Hi contains Hubbard operators that acquire a non-zero mean value, our lattice 
model is mapped into an impurity model with an external field determined self-consistently by equation (|29p . 



9 



Our final step consists in finding a self-consistent equation that determines the hybridization A(z) in terms of the 
lattice Green's function. Comparing Eq. (|19|l for Gy (after making My = i,jM ) with the equivalent expression for 



(0) 



cg } (tw n ) = %M (jw n ) + 2 MoftWnjEttGg 5 («*,) for i, j ± 0, 

fe#0 



(30) 



we obtain for i = j = 0: 



(iw n ) 

k,l 

1 - 53 E ofc GS ) (^ n )E io M («w„) 



= Mn 



k.l 



(31) 



Using equation (|28[1 for A.(iui n ) and the Fourier transform for the local Green's function, we obtain the self-consistent 
condition that determines the hybridization, 



1 53 [M- V„) - E(k)]- 1 = [M- V«) - A(iu; n )] 1 . 



(32) 



Formally, Eq. Ij32(l closely resembles the standard form of the DMFT self-consistency condition!. We should, 
however, keep in mind the fact that this is a matrix equation involving Green's functions for the Hubbard operators. 
To solve a problem described by the Hamiltonian (lll.'il) within this scheme, implies to solve the impurity problem, 
B24J1 . find an irreducible cumulant matrix Mo, then use the self-consistency condition (|32(l together with equation 
(|29|1 to determine the new hybridization, A(iw n ), and the new field hi mp . The process is repeated until convergence 
is reached. 



A. The Hubbard Model 



In order to get a feeling of how this approach works, as well as to check its validity, we can do a simple exercise by 
applying the procedure to the Hubbard model, for which we should be able to recover the known DMFT results. In 
this particular case, the Hamiltonian will be 

Hhub =5>M*; T +^ U ) + (^-2mR 22 ]+EE UjiXfX^ + X^Xp + aX^Xp+aX^X^}, (33) 

i i,j <7 

where \i represents the chemical potential, U the on-site Coulomb interaction and ty the hopping matrix elements, 
which are equal to —t if i and j are nearest neighbors and zero otherwise. We can introduce the hopping matrix, 



E(k) = I E 2S °"(k) E 2S ff2 (k) ] = [ ae k e k 



(34) 



where e k is the Fourier transform of iy . Next, we introduce the Green's function matrix and the irreducible cumulant 
matrix: 



G 



£^0cr crO £^0cr 2(7 
£^(x2 crO 2(T 



M 0<T <t0 M 0a 2<T 
M = I M s2 a0 M s2 2S 



(35) 



Our goal is to show that the matrix equation i|32|) reduces in this particular case to the well known DMFT self- 
consistency condition, 



G l (iu) n ) = iuj n + fx - A Q (iu} n ) 



1 E 



N 4^ iuj n + fx - e fe - £(ia>n) 



(36) 
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where E(iw„) represents the self-energy of the single-particle electron Green's function and Ao(iw„) the hybridization 
of the effective single-impurity problem. The electron single-particle Green's function, Q, is given by a sum of all the 
elements of G, as shown in Eq. 122(1 . Similarly, we define the quantity 

M{iu n ) = M 0a a0 (iu n ) + M s2 2s {iuj n ) + a[M 0a 25 {iu n ) + M s2 M {iuu n )]. (37) 

ft is straightforward to show that the matrix Green's function has the structure: 

G(k, iuo n ) = [M- 1 ^) - E(k)]- 1 = L_— M(ta; n ) + 5- £fc (a -1 ) , (38) 

where 8 = det{M}. Using Eq. ((38(1 and the expression ((22(1 for the electron Green's function, together with the 
definition of M, we have 

g(k, iu n ) = 1 , (39) 

M L (iLLi n ) - e k 

which shows that M.(iu) n ) is the single particle irreducible cumulant, .A/f _1 (iw„) = iuj n + fi — E(ta; n ), where £(iw ra ) 
is the single particle self-energy. On the other hand, making the summation over k in Eq. ((38(1 and introducing the 
result in the self-consistency relation ((32() gives: 

la 

A = A I a 1 I , (40) 



with 



b 1 1 , *(o-l) 
A = — , a = - > , b=— -. (41) 

k 

Taking into account the structure of A given by Eq. 1(40(1 . the impurity problem described by the statistical operator 
((2411 reduces to an Anderson single impurity problem with a hybridization A . Finally, we have 



Ao(i<jj n ) = M 1 (iuj n ) 
which represents the standard DMFT self-consistency equation 



N ^ M~ l {iu n )-6 k 



(42) 



IV. A FUNCTIONAL APPROACH 

Functionals of Green's functions provide an elegant and powerful tool for formulating dynamical mean-field the- 
ories. Within this framework, the DMFT equations are obtained as a direct result of a simple approximation on 
the Baym-Kadanoff functional 1 . The extended dynamical mean- field equations have also been rigorously derived 
using a functional technique^. Moreover, the functional perspective seems to offer the natural framework for cluster 
generalizations of the dynamical mean-field theory that include effects of short-range correlations. The purpose of 
this section is to introduce a strong-coupling generalization of the functional DMFT formulation- . The construction 
of a generalized Baym-Kadanoff-type functional will allow us to present an alternative derivation of the generalized 
DMFT equations, and, more importantly, will set the basis for a future cluster extension of the G-DMFT approach. 

The main idea of the functional interpretation of DMFT is to formulate a local approximation for the Baym- 
Kadanoff functional $[G], which is the sum of all vacuum-to- vacuum skeleton (two-particle irreducible) diagrams 
constructed with the full propagator G and the interaction vertices. This functional has the property that 



5$ 

5Gji{iw r , 



T,ij{iu) n ) = ^ — , , (43) 



where Tnj{iui n ) is the self-energy. As the coordination number of the lattice goes to infinity, the self-energy becomes 
local and $ will depend only on the local Green's functions, Gu- 



$ = 53 c/>[G u ], (44) 
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where 4>[Gu] is a functional of the local Green's function at site i only. This scheme cannot be transferee! directly to 
approach our strong-coupling problem. As we mentioned above, within our strong-coupling perspective, the natural 
quantity to describe the system is the two-point irreducible cumulant, Mij(iui n ), rather than the self-energy. Our task 
is to identify a conjugate quantity (which cannot be the Green's, function but rather something with dimensions of 
energy), and to construct a functional of this quantity that will determine M by a relation similar to l|43[) . 

Following Refs. and^, we start by introducing two time dependent external sources, one , Jt/( T ~ T ')i coupled 

to the X-operator Green's functions, and the other, A"^(r), coupled to the average = (X" 13 ) of the bosonic 
Hubbard operators that occur in the non-local Hamiltonian H±. The corresponding generating functional will be 



Q[J,X] = --lnTrp[J,X], 

P 



(45) 



where the statistical operator p is given by 
p[J,X] = e- pHa f exp 



(46) 



The field J-^ (t — r') is equivalent with introducing frequency dependent coupling constants E^{iuj n ) 



J* 



Jj^(iu! n ). The Green's functions can be expressed as functional derivatives of the generating functional, fi[J, A], with 



respect to the frequency dependent field, 



Gg(io; n )-f3(X?)(X?)=-l3 



8Sl[J, A] 



(47) 



J=0,A=0 



where the average value Qf = (X?) is non-zero only if X? is a Q-operator. The generating functional fl[J, A] reduces 
to the grand-canonical potential, il, in the absence of the external sources. 

The next step is to integrate Eq. (|47|l . As we are interested in a strong-coupling equivalent of the Baym-Kadanoff 
functional, we will express the Green's function in terms of the two-point irreducible cumulant, rather than the 
self-energyi£. Using the notation A for a tensor with elements A^ , we can formally write equation 1)19(1 as 



G= I ME M, 

where I is the identity tensor, 1^ = SijS^v By integrating (|4*7ll with G given by we obtain 

Q[J, A] = (l + i Tr + Tr[M£] + E[J, A], 



(48) 



(49) 



with Slo the grand-canonical potential in the atomic limit. The quantity £ is the renormalized coupling constant (or 
renormalized generalized hopping) and can be represented diagrammatically as a sum of chains of "hopping" lines 
and irreducible cumulants. Analytically we have 



p,q a,/5 

The last term in Eq. (|49|l is a functional with the property 



E(I — ME)" 1 



(50) 



Tr 



M- 



S£ 



(51) 



Also, because the functional derivative of f2[J, A] with respect to A gives the average value of Af , we will have 

si 



SE 1 
= Tr 



M- 



5X^ 



(52) 



The standard way to proceed is to perform a Legendre transform of the generating functional with respect to the 
interacting single-particle Green's function and to eliminate the external source J in favor of G in the functionals. 



12 




FIG. 4: Generalized skeleton diagrams contributing to the functional H*. The contribution of the external field A is represented 
only in diagram (i) by dotted lines. The numbers in parentheses represent the symmetry factors (if different from 1) of the 
corresponding diagrams (the number of ways in which a graph can be made isomorphic with itself) . The renormalized coupling 
£ is represented by double lines. 



Instead, by analyzing equations Ij51|l and 1)52 jl. we identify the renormalized coupling constant £^ as the quantity 
conjugate with the irreducible cumulant and use it as an independent variable. Consequently, we introduce the new 
functional fl* [£, A; E] and S* [£, A; E] , with 

Cl*[£[J, X], X; E + J] = fi[J, A], H* [£[J, \],\;E + J] = E[J, A] . (53) 

The dependence on the last variable, E, involves only the components of the coupling constant tensor corresponding 
to Q-operator indices and we have 



SE* 
6E» v 



m<Xj)- (54) 



The functional derivative with respect to A of both SI* and S* will give the average value of the bosonic X operator, 
while for the Luttinger- Ward-type functional S* we have 

5£ j f{iuj n ) 

This equation is the strong-coupling equivalent of l|43|) and determines the irreducible two-point cumulant as a func- 
tional of the renormalized "hopping" matrix, the bare diagonal coupling elements and the external fields. The 
stationarity condition SQ*/S£j- t — yields the relation between the renormalized couplings and the cumulants, 

£%y = [E(I — ME) -1 ]^", which is identical to the definition Q50[l. Again, at stationarity and in the absence of the 
external fields the functional f2* reduces to the grand canonical potential. 

For the case with no Q-operators in the non-local Hamiltonian our construction is complete. In this particular case 
the only relevant dependence in the functionals is on the renormalized coupling £ . The functional f2*[£] is formally 
analogous to the functional r„ eu) [G] introduced by Chitra and Kotliar— from a weak-coupling perspective. In our 
picture the strong coupling correspondent of the self-energy is, as noted before, the irreducible cumulant, while the 
renormalized and bare coupling constants correspond to the full and non-interacting Green's function, respectively. 
Also, the strong-coupling analog of the Luttinger- Ward functional $[G] is E*[£}. Similarly to its weak-coupling 
counterpart, E*[£] can be expressed as a sum of all vacuum-to-vacuum skeleton diagrams containing bare cumulants 
and renormalized "hopping" lines. All the skeleton diagrams are two-particle irreducible, i.e. they cannot be separated 
into disconnected parts by cutting two £ lines. 

In general, however, the diagrammatic expression of EE* is more complicated and involves, in addition to the bare 
cumulants, both renormalized coupling constants, J/, and bare coupling constants, . Also, lines corresponding 
to the external field A have to be added. The generalized skeleton diagrams have to satisfy two conditions: i) Each 
diagram is two-particle irreducible with respect to the renormalized "hopping"; ii) By cutting any bare "hopping" 
line the diagram is divided into two disconnected parts. These conditions determine a tree-like structure for the 
generalized skeleton diagrams with skeleton blocks coupled by bare E lines. This structure is illustrated in Fig. 0] 
for a few low-order contributions to EE* . Due to the tree-like structure of the diagrams, it is rather difficult to find 
a local approximation for the functional H*. In order to overcome this problem, we proceed to the last step of our 
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construction and define the Legendre transform of fl* with respect to the average of bosonic Hubbard operators, 

T[£,Q]=n*[S,\;E]-Y / E ( 56 ) 

i 1-ieBq 

where the bare coupling E?? is treated as a parameter and is not written explicitly as an argument of the functional 
r. Using the condition — = 5T/6Q^ to eliminate A, we obtain 

Y[£,Q] = n + -Y,Tr Hi-m[£MM + \Y, Tr ^- e ^-\Y. E QtE^Q^ + n^Q]- (57) 

The properties with respect to £ of the new functionals are similar to those of f2* and 3*, namely the stationarity 
condition ST/SS^f 1 = yields equation l|5U|) . and we have 



In addition, stationarity with respect to Q, 5T/6Q^ — 0, yields 

The diagrammatic representation of the new Luttinger- Ward-type functional '5 contains skeleton diagrams that are 
two-particle irreducible with respect to the £ lines and have cumulants dressed with EQ contributions. The structure 
of the diagrams is illustrated in Fig [3] Using the diagrammatic representation of the functionals $ and 3*, as well as 
that of Q, and taking into account the symmetry factors associated with each diagram, it is straightforward to check 
that 

*[£,Q] = -*[£, A = 0;£] + iE E Qf^Qi' ( 60 ) 

i,j /i.weBo 

where Qf = (Xf ) = 5E*/S\? | A=0 . 

We are prepared now to introduce the generalized DMFT from a functional perspective. By analogy with the 
standard approach^, we construct an approximate theory by making the ansatz that the Luttinger- Ward functional 
\& depends only on local renormalizcd "hopping" matrices, i.e., 

*[£>Q] = 52l>[£u,Q], (6i) 

i 

where ip is a functional that depends only on the local renormalized coupling constants at site i. In diagrammatic 
terms, it means that all the bare cumulants contained in a skeleton diagram contributing to tp are defined on the same 
site. A direct consequence of this ansatz is that the two-point irreducible cumulant of the theory is local, 



5* 

8£ il {iu) n ) 



Mij(iuj n ) = — - Sij, (62) 



with Sij being the Kronecker S symbol. The stationarity condition for T yields 



f — 



E(i - ME)- 1 ] = i E E W( r ^ ME ( fc )) _1 > ( 63 ) 



N 

k 



where M is given by (|62|l and, for the second equality, we assumed that the system is homogeneous. Relation 16311 
represents our generalized DMFT equation and it should be viewed as a functional equation for the local renormalized 
coupling matrix. To cast it in a more familiar form, let us observe that using the expressions for the Green's function, 
Eq. I|48|). and for the renormalized coupling, Eq. (JSDJ, we have 



G = M£M + M. 



(64) 
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FIG. 5: Example of diagrams contributing to ^[£,Q\. The average of a diagonal operator, Q™, is represented by a square. The 
symmetry factors are written in parentheses. If = 0, the remaining basic skeleton diagrams (a, e, g, ...) are identical with 
the remaining diagrams in Fig. |H(a, b, c, ...), and ^ = S*. In general, we have the basic diagrams plus diagrams obtained by 
decorating the bare cumulants with lines and squares representing E^YQ^. 

Taking into account the local character of the irreducible cumulant, we obtain, for a homogeneous system, 

-l 

(65) 

where Gu{iuj n ) has to be understood as a functional of £. 

To be able to solve Eq. (|63fl or Eq. I|65|> we need to know the functional dependence of ^, at least approximately. 
A direct approach is rather hard in practice, therefore it is more convenient to generate the functionals from a purely 
local theory. A simple inspection of the diagrammatic expression for <3> (see Figure |SJ, shows that, in the case when 
all the bare cumulants are defined on the same site i — 0, the same diagrams will correspond to an impurity problem 
defined by the statistical operator given by Eq. 124f> . In fact, taking into account the previous construction of the 
strong-coupling functionals, we have 

T imp [D, Q] = ^r + \^Tr ln[I - M imp [V, Q] A] + ~ E Tr[M imp [V, Q]V] (66) 

where T> = A [I — Mi„ lp A] _1 is the renormalized hybridization, = (X^) and A Ml/ (0) is the hybridization at zero 
frequency. In addition, we have to keep in mind that the bare impurity cumulants are dressed by external field lines, 
ftf . The diagrammatic expressions for ^>[£ , Q] and ^i mp ['D, Q] are identical and we have 9/N — ^>i mp provided 

= J2 E ZjQi = h im P + E ^"(0)Q". (67) 

3," v 

Equations (|67|l represent the self-consistency condition for our strong-coupling generalized dynamical mean-field the- 
ory. The first of these equations can be rewritten in the more familiar form, Eq. Il.'j2[| . by simply using the expressions 
of the renormalized quantities £ and T>. A direct consequence of the first self-consistency condition is the equality of 
the irreducible cumulants, M^(iw„) = Mj mp (io; n ). Finally, let us mention that this correspondence with the impurity 



G H {iw n ) = i E 
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problem is also useful in establishing a relation between the grand-canonical potential of the lattice and the grand- 
canonical potential of the impurity. At stationarity both T imp and T reduce to the corresponding grand-canonical 
potential, while *ff/N = ^i m p- Therefore we have 



N 



1 \- 



h3 



Tr\n[I — MA] — ^ Tr ln[I — ME(fc)] 



(68) 
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APPENDIX A: ALTERNATIVE DERIVATION OF THE G-DMFT EQUATIONS 

The purpose of this Appendix is twofold: 1) To present an alternative formulation of our strong-coupling theory 
that is closely related to the more familiar perturbation expansions based on Wick's theorem and translate the basic 
quantities involved in this approach in the corresponding language. 2) To propose a scheme that allows one to 
expand around the DMFT solution. The central idea is to decouple the nonlocal term using a generalized Hubbard- 
Stratonovich transformation and to trade in the X-operators for a set of canonical fcrmionic and bosonic fields. In 
terms of these new auxiliary fields, the standard techniques of the perturbation theory can be applied. This formulation 
represents a generalization of the strong-coupling expansion used by Pairault, Senechal, and Tremblayi^ for the study 
of the Hubbard model. 

Starting with the partition function 



Z = Tr{ e pH ° f exp 



fdr ]T Y, X t(T)EgX!(T) 

J ° <i,j> \i,v 



(Al) 



we perform a Hubbard-Stratonovich transformation and express the non-local part as a Gaussian integral over the 
auxiliary field y , 

Z = Z [V( exp J f dr £ £ Cf (r)(Xo W " S ™® \ > ^ 
J {■' <ij>w J 

with 

Si„t[C] = -]»(7>M><{ / <h-y (-r)AT(r; Y,"(rK,'Vi J> ) ■ (A:!, 



\n(f exp|-^dr^[Cf(r)Af(r)+Af(r)Cr(r)]| ^ . 



The field contains both fermionic (Grassmann variables) and bosonic (complex) components, as determined by 
the corresponding X-operators. The free propagator appearing in the quadratic term of Eq. (|A2f> is Xo% = ®^ . 
Notice that the auxiliary fields have canonical fermionic or bosonic statistics and all the complications related to the 
unconventional algebra of the Hubbard operators are buried in the average (. . .)o with respect to the unperturbed 
ensemble. Due to the nature of the local Hamiltonian the interaction term Si n t[£] can be written as a sum of 
purely local contributions, each having an expression similar to Eq. (|A3(I but without the summation over the site 
index i. Expanding in the C-fields, we have explicitly 

oo 

*WC]=EE^ E C7f( Ml ,..,Mn;n,...,r„)Cf™(r n )...cr(n), (A4) 

where C" is the n-point bare cumulant defined by Eq. (|14f> . We can now develop a standard perturbation theory 
for the C-neld, having xo = E as the free propagator, and vertices given by C™. The relation between the Green's 
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function of the X-operators and that of the auxiliary fields can be easily obtained if we notice that, due to the form 
of the interaction term (|A3() , we have 

G^,j;t,t') = iy©C $00^ exp{-So[C]}, (A5) 
where Sq[Q represents the quadratic term in Eq. I|A2|I . After performing two integrations by part, we obtain 

G = -Xo 1 +Xo 1 XXo\ (A6) 
where we used the tensor notation, and we introduced the full propagator, \i f° r the auxiliary fields, defined as 

X?/(t,t/) = -(f + (CfXCj), (A7) 

where T is the (imaginary) time ordering operator and (...) represents the average over the statistical ensemble. In 
terms of the self-energy for the C variables, £^ = Xo 1 — X^ 1 > onc nas 

G = (S c - 1 -Xo)^ 1 , (A8) 

which shows that corresponds to the irreducible cumulant M of the X-operator formulation. Also, the full Green's 
function x = (x^ 1 — ^c) _1 = ^(1 — ME) -1 is the renormalized coupling costant £ of the original theory. In addition, 
in the presence of bosonic Hubbard operators with non-zero mean value, we have 

(cr>=E B £W>> (A9) 

showing that the average of the auxiliary field represents the effective "magnetic" field of the original formulation. 
This equation completes the correspondence between the the formulation in terms of X-operators and that in terms 
of auxiliary ^-fields. 

It is now straightforward to re-derive our generalized mean-field equations using the standard techniques of a 
perturbation theory for the auxiliary fields. For example, the Baym-Kadanoff functional, $[x], can be defined in the 
usual way as the sum of all vacuum-to-vacuum skeleton diagrams. The corresponding free-energy functional can be 
written as 

A Ix, Ql = ^E Trln fa 1 - [x, Q]] + \ E Tr $dx, Q]x\ \ E E Qi (xo)f/QJ + $bc, Q], (Aio) 



where the possibility of a non-zero average, Q, for the X-operators has been explicitly considered. Taking into account 
the correspondence between the auxiliary field formulation and the original X-operator formulation, we can see that 
A[v, Q] is the equivalent of the functional T[£, Q], while the Baym-Kadanoff functional <5[x, Q] corresponds to &[£ , Q], 
up to an £- independent term equal to the atomic grand-canonical potential Qq. 

We can construct now a strong-coupling expansion to a certain order in (Eij/U), as in the work of Pairault, Senechal, 
and TremblayJ- . However, any finite order expansion will intrinsically contain the problem of the exponentially large 
degeneracy of the atomic ground state, leading to the break-down of the solution in the low temperature limiljiS. 
Next, we can recover our strong coupling DMFT scheme as a local approximation for the Baym-Kadanoff functional. 
This renormalized perturbation approach avoids the problem of the low temperature barrier. Further, we can exploit 
the idea of introducing auxiliary degrees of freedom and construct a formulation which, in principle, allows one to 
determine corrections to the DMFT approximation. The basic idea is to introduce the auxiliary fields in such a way 
as to obtain the DMFT result as the lowest order approximation. Starting again from the partition function IjAljl . 
instead of performing directly a Hubbard-Stratonovich transformation, we first add and subtract in the exponent a 
term equal to 

f dr f dr'^^Af(r)A^(r-r')Ar(r'), (All) 

where A^ ly (r — r') is the hybridization function for the impurity problem associated with the DMFT approximation. 
Notice that this is a purely local quantity. Next, we introduce the auxiliary field £ that decouples the quadratic term 
containing xo( T ~ r ') = ES(t — t') — A(t — r'). The effective theory for the new variables will have a free term 



«= - I dr [ dr' E E^HXo 1 )^-^)^'), (A12) 
• /:: 7 ' <i,j>w 
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and an interaction term 



S int [£] = - In (f cxp j - jf " drdr' XfAgX? + &X? + X?& j ^ , (A13) 

where, for simplicity, we omitted the summations and the imaginary time dependence. Again, the interaction part 
can be written as a sum of purely local terms. However, the vertices of the new theory are no longer the bare 
cumulants of the Hubbard operators. Instead, let us notice that each term Si„t[£,i] of the interaction part represents 
the generating functional — ln[Tr({/3j mp }] for the impurity problem defined by Eq. (|24|) . In the present formulation, 
the external fields h°^ p are represented by the average values (£ Q/3 }. A direct consequence of this correspondence is 
that the second order cumulant generated by Si n t is equal to the impurity Green's function Gi mp associated with 
the generalized DMFT solution of our problem. Further, the relation between the propagator for the £-field and the 
X-operator Green's function can be derived as before and we have 

G=(t^-Xo)-\ (A14) 

with Ej the self-energy for the £-field. 

Let us construct now the lowest order approximations for the pertubation expansions formulated in terms of 
auxiliaryy fields C an d £, respectively. In the first case, the simplest diagram contributing to E<; is equal to the 
atomic Green's function. The corresponding result represents the Hubbard-I approximation 19 . In the second case, the 
zero-order approximation for the self-energy is given by the DMFT impurity Green's function, E^ = Gi mp . Using 

the correspondence relation (| A14|l with xo — E — A, as well as the equation G^ mp {iui n ) — iu>„ — A(iuj n ) — Ej mj) (m; n ), 
we have 

G(tw n ) = (G^ p (iu n ) - E + A)- 1 = [iu n I - E Jl imv {iuj n )}- X . (A15) 

This shows that, in the formulation using auxiliary £ fields, the generalized DMFT result is recoverd as the lowest 
order approximation. Obviously, in order to determine the hybridization function A, one has to first find the solution 
of the problem within the DMFT scheme. However, the strong coupling expansion constructed here allows us to 
determine in a systematic way the corrections to the DMFT result by computing higher order contributions to the 
self-energy Eg. 
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